Extensibility enables locomotion under isotropic drag 
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Anisotropic viscous drag is usually believed to be a requirement for the low Reynolds number 
locomotion of slender bodies such as flagella and cilia. Here we show that locomotion under isotropic 
drag is possible for extensible slender bodies. After general considerations, a two-ring swimmer and a 
model dinoflagellate flagellum are studied analytically to illustrate how extensibility can be exploited 
for self-propulsion without drag anisotropy. This new degree of freedom could be useful for some 
complex swimmer geometries and locomotion in complex fluid environments where drag anisotropy 
is weak or even absent. 



Due to the absence of inertial forces, low- Reynolds 
number locomotion is subject to interesting mathemat- 
ical and physical constraints [1-3]. In particular, loco- 
motion by time-reversible strokes is ruled out by Pur- 
cell's scallop theorem [4, 5]. To escape these constraints, 
microorganisms swim by either propagating deformation 
waves along slender appendages, termed flagella, or ro- 
tating them. Anisotropic viscous drag is believed to be 
the fundamental property enabling drag-based propul- 
sion of slender filaments [3, 6, 7]. It is a classical result 
that for a slender filament moving in an unbounded New- 
tonian fluid, the Stokes drag is almost twice when mov- 
ing perpendicular than parallel to its axis [3, 8]. This 
drag anisotropy allows propulsive forces to be created 
perpendicularly to the deformation of the filament. Un- 
der isotropic drag, it is generally accepted that locomo- 
tion of this kind would be impossible [6, 7, 9-11]. 

Unlike in Newtonian flows, drag laws in more complex 
media, and their consequences on locomotion, remain 
largely unexplored. Theoretical studies, via Brinkman 
models, suggest that porosity enhances drag anisotropy 
[12, 13], explaining, e.g., the increase in propulsion speed 
of C. elegans in a granular medium. Recent experiments 
also measured and characterized granular drag in beds of 
glass beads and granular media [14-17], which have been 
applied to study locomotion in sand [18]. Besides the 
fluid medium, the geometry of a swimming body also 
plays a role in the drag law. Some flagella, such as 
those of Ochromonas, possess rigid projections termed 
mastigonemes, protruding into the fluid [1]. In these ge- 
ometries the viscous drag in the longitudinal direction of 
the flagellum is increased, possibly resulting in a more 
isotropic drag. In situations where drag anisotropy is 
weak or even absent, what are the alternative mecha- 
nisms, if any, offered by physics to achieve locomotion? 
In this letter, we point out a new degree of freedom en- 
abling inert ialess swimming, namely extensibility. Using 
a general derivation and two simple geometrical models, 
we demonstrate that the periodic stretching and con- 
traction of a filament allow self-propulsion even under 
isotropic drag. 



We start by considering the general calculation of 
Becker et at. [7] showing that drag anisotropy is required 
for the propulsion of inextensible swimmers. Here we re- 
visit their derivation by relaxing the inextensibility con- 
dition. Consider a filamentous swimmer of total length 
and denote by r(s,t) the instantaneous position of 
material points along the filament, where s is the ar- 
clength. The time rate of change of the average swimmer 
position, f (t), is given by 
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where we have denoted L = dL/dt. Under isotropic drag, 
we have dr/dt = u oc f, where f is the local drag force 
per unit length, and thus dr/dtds oc ids = for 
force-free swimming. We therefore find that the second 
term in Eq. (2) disappears, and thus 
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For an inextensible swimmer, L{t) = 0, leading to 
dr/dt = 0: no net propulsion is possible under isotropic 
drag [7]. As a difference, for an extensible swimmer 
{L 7^ 0) no general conclusion can be drawn, suggesting 
the relaxation of the drag anisotropy requirement when 
extensibility is permitted. We demonstrate below how 
this additional degree of freedom can be exploited for 
self-propulsion by considering two simple examples. 

As the first example, consider the motion of two ex- 
tensible slender rings (red solid and blue dashed circles 
in Fig. 1) connected by two rotating rigid rods (each of 
length 2r). As a result of the rod rotation, each ring 
expands and contracts periodically, tracing in time the 
surface of a torus of inner radius R and a circular cross- 
section of radius r (see Fig. la). In Cartesian coordi- 
nates, the motion of each ring r{t) = [x(t), y{t)^ z{t)] can 
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FIG. 1: (a) Notation for the extensible two-ring swimmer, (b) Schematic illustration of its cyclic deformation over one period, 
T (see text for details). 



be described as 
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where, G [0,27r], uj = 27r/T is the angular frequency, 
T is the period of the motion, and (j)o is the phase of 
the motion. We non-dimensionalize lengths by the in- 
ner radius of the torus time by l/cj, and the di- 
mensionless deformation kinematics of the two rings, 
fi,2(?) = [^i,2W,^i,2(?),^i,2(^)], are given by 

xi,2{i) = {1 + f [1 - cos{i + (^1,2)] } cos6>, (7) 
yiAt) = {1 + r [1 - cos(? + (/)i,2)] } sin^, (8) 
^1,2^ = rsin(t + (^i,2), (9) 

where we assign ^1 = (red solid ring) and 02 = tt 
(blue dashed ring). All tilde variables in this letter are 
dimensionless. 

We illustrate the cyclic deformation of this swimmer 
over one period in Fig. lb. Hydrodynamically, material is 
being created/destroyed when a ring extends/contracts, 
introducing an additional mechanism for varying the hy- 
drodynamic drag. The non-zero translational velocity 
expected to arise is a direct result of the imbalance of 
vertical viscous force due to the difference in the total 
arc-length of the two rings. Since the vertical motion 
and variation of the circumference of the two rings are 
out-of-phase, the vertical velocity of the ring with a larger 
circumference always points in the same direction (in the 
case considered here, in the —z direction). Therefore, 
we expect a net unidirectional vertical force, and hence 
swimming, in this direction. Note that when the two 
rings have, instantaneously, exactly the same diameters 



(at i = (2n + l)7r/2, where n = 0, 1, 2, ... ; see the up- 
per right and lower left panels in Fig. lb), the viscous 
forces acting on the rings balance, and no swimming is 
produced. Note also that if there was only one ring, 
say fi, the vertical motion of the ring, would be exactly 
canceled by translational swimming velocity at the same 
speed, resulting in no apparent motion. 

Quantitatively, we apply a general local drag theory 
to this two-ring swimmer. We neglect the hydrodynamic 
impact of the rods connecting the rings. The local vis- 
cous force density acting on the filament is then given by 
f = — [^||tt + ^^(1 — tt)] • u, where t is the local tangent 
vector and u is the local velocity of the filament. The 
drag coefficients and characterize the motion of a 
slender rod parallel and perpendicular to its axis respec- 
tively. Their specific values depend on the geometry of 
the rod and properties of the fluid medium. For a slightly 
distorted slender filament in an unbounded Newtonian 
fluid. Gray and Hancock [6, 9] derived explicit analytical 
expressions for the drag coefficients, which were later im- 
proved by Lighthill [19]. For more complex environments, 
these drag laws remain mostly unknown. Under isotropic 
drag, we have ^± = [20-23]. Here, we keep their val- 
ues general in the calculations, and show that the final 
swimming speed of the two-ring swimmer is independent 
of these coefficients. The swimmer hence works equally 
well under any local drag law, including isotropic drag. 

The local velocity along a ring, u = dr/dt-\-\J -\-ftxr^ 
is composed of two parts: the prescribed deformation 
velocity, dr/dt, and the unknown swimming velocities 
U + X r to be determined. At low Reynolds num- 
bers, the total force and total torque on a swimmer 
have to vanish. Hence, we have X^i^i lo'^^"^ ^idsi = 
Xli=i Jo'^^"^ 'Tidsi = 0, where r = rxf is the local viscous 
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FIG. 2: Time- variation of the dimensionless swimming veloc- 
ity of the extensible two-ring swimmer, Uz, for different values 
of r. 



torque density, leading to six equations to determine the 
unknown swimming kinematics (U,r2). By symmetry, 
we have = 0, and Uz is the only non-zero translational 
velocity component. After some algebra we obtain the 
dimensionless swimming velocity 
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We display the variation of the swimming velocity as 
a function of time for different values of r in Fig. 2. 
We obtain unidirectional swimming {Uz > 0) at a time- 
averaged speed independent of the value of the drag co- 
efficients, (Uz) = 7^^/2(1 + f): this extensible two-ring 
swimmer can therefore self-propel under isotropic drag. 

This idea can be extended to more complicated geome- 
tries consisting of a curved structure built upon another 
curved structure, for instance, toroidal helices (a helix 
built upon a circle) and superhelices (a helix built upon 
another helix). A toroidal helix has been studied as an 
idealized geometrical model for dinoflagellates [24, 25]. A 
dinoflagellate can be propelled by propagating a toroidal 
helical wave, where the kinematics of such a wave im- 
plicitly assumes extensibility due to the intrinsic length 
differences in the geometry. Consider a toroidal helical 
wave with amplitude r, wave- number k = 27r/A and an- 
gular frequency co', propagating along a circle of radius R 
(see Fig. 3 inset). When lengths are non-dimensionlized 
by and times by l/cj, the dimensionless kinematics 
R = [x' ^y' ^ z'] can be expressed as 

x'{s',t) = Rcos{s'/R)^rcos{s' -i)cos{s'/R), (11) 
y{s,i) = Rsm{s/R)^rcos{s-i)sm{s/R), (12) 
z\s\i) = ±rsm{s' -i), (13) 

where s' G [0, 27tR] parametrizes the toroidal helix, and 
where all tilde parameters are dimensionless. Note that 
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FIG. 3: Dimensionless translational velocity Uz as a function 
of the dimensionless toroidal wave amplitude, f = r/c, under 
isotropic = ^^^^ solid line) and anisotropic = 2^||, 
red dashed line) drag, for the case of = 8. Inset: notation 
for a toroidal helix swimmer (see text for details). 



the dimensionless radius R takes only integer values for a 
closed toroidal helix. The ± sign represents different chi- 
rality of the helix. We apply the same local drag model as 
above. In order to make analytical progress, we consider 
the small-amplitude limit f <C 1 and perform asymptotic 
expansions in powers of f. Swimming occurs at O(f^) 
and takes the dimensionless form 
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While the rotational velocity ft vanishes under 
isotropic drag {(,±/(,\\ = 1), the translational velocity U 
along the torus axis is independent of the drag coeffi- 
cients. This is due to the intrinsic difference in length 
(and hence drag) embedded in the curved geometry of 
toroidal helices. The toroidal helix has shorter lengths 
on the sides closer to the center of the circle, creating an 
overall imbalance of hydrodynamic drag. It is interesting 
to compare two physical limits, namely when there is no 
drag anisotropy = ^||) when there is maximum 
drag anisotropy for an asymptotically slender filament in 
a Newtonian flow = 2£^\\)- Using these values in the 
general local drag theory, we find that the axial propul- 
sion velocity in these two limits are indistinguishable for 
small f (Fig. 3), as predicted by the asymptotic analysis. 
As r increases, the drag anisotropy increasingly enhances 
propulsion (up to ~ 30% when f = 5) but extensibility 
alone, under isotropic drag, still enables swimming with 
significant magnitudes. A similar example that could ex- 
ploit extensibility to produce swimming under isotropic 
drag are super helical waves, which serve as a geometrical 
model for insect spermatozoa [26]. 
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In conclusion, in this letter we have revisited the re- 
quirement of drag anisotropy to achieve self-propulsion at 
zero Reynolds number. We demonstrated expHcitly, via 
two simple swimmers, that extensibility provides a mech- 
anism for swimming under isotropic drag, which might 
be relevant for some complex swimmer geometries and 
motion in porous or other complex media where drag 
anisotropy is weak or absent. As a practical side- note, 
the two-ring swimmer described above could be actuated 
experimentally by external rotating magnetic fields if fer- 
romagnetic materials, such as metal nanowires [27, 28], 



were used for the phantom rods. 
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